CMB Neutrino Mass Bounds and Reionization 
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Current cosmic microwave background (CMB) bounds on the sum of the neutrino masses assume 
a sudden reionization scenario described by a single parameter that determines the onset of reion- 
ization. We investigate the bounds on the neutrino mass in a more general reionization scenario 
based on a principal component approach. We found the constraint on the sum of the neutrino 
masses from CMB data can be relaxed by a ~40% in a generalized reionization scenario. Moreover, 
the amplitude of the r.m.s. mass fluctuations erg is also considerably lower providing a better con- 
sistency with the low amplitude of the Sunyaev-Zel'dovich signal recently found by the South Pole 
Telescope. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 



INTRODUCTION 

The high precision measurements of Cosmic Microwave 
Background (hereafter, CMB) anisotropies made by the 
Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite have provided not only a wonderful confirmation of 
the standard model of cosmological structure formation 
but also relevant information on key parameters in par- 
ticle physics. One example is the sum of the neutrino 
masses and the neutrino mass hierarchy. 

The most recent data release from WMAP after seven 
years of observations presented a bound on the total neu- 
trino mass of Y>m v < 1.3cV at the 95% c.l. 1]. This 
bound is approximately a factor five better than the cur- 
rent laboratory experimental upper limit inferred from a 
combination of beta-decay experiments and neutrino os- 
cillation data (see e.g. 0)- The CMB bound on neutrino 
masses is also considered the most conservative limit from 
cosmology. Indeed, including information from galaxy 
clustering and luminosity distance data the constraint 
can be further improved to Sm„ < 0.55eV at 95% c.l. 

while a limit of Y>m v < 0.28eV at 95% c.l. can be 
obtained by including redshift-dependent halo bias-mass 
relations [3]. 

It is however important to be aware of the theoretical 
modelling behind the constraint based on cosmological 
measurements. A model of structure formation based on 
dark matter, adiabatic primordial fluctuations and dark 
energy is assumed and the removal of one of these as- 
sumptions can in principle affect the CMB limit. For 
example, the inclusion of isocurvature perturbations 
dark energy [l[ or modified gravity Q can all relax the 
CMB upper limit on the neutrino masses. 

In this Brief Report we investigate another possible 
theoretical caveat that could affect the CMB bound on 
the sum of the neutrino masses, i.e. the modelling of the 
reionization epoch. It is often assumed in the current 
cosmological data analysis that reionization is a sudden 



event at redshift z — z re , i.e. this process is usually 
described by a single parameter with the free electron 
fraction x e increasing from ~ 10~ 4 up to 1 for redshifts 
z < z re (~ 1.08 for z < 3 when taking into account 
Helium reionization). While this scenario can properly 
describe several reionization scenarios, it can't obviously 
describe more complex reionization scenarios as for ex- 
ample double or not-monotone reionization. Given our 
current ignorance about the thermal history of the uni- 
verse at redshifts z > 6 it is important to consider all the 
possible reionization scenarios allowed by data when de- 
riving the most conservative constraint on a cosmological 
parameter such as the sum of the neutrino masses. 

Here, we indeed assume a more general reionization 
model following the principal components method sug- 
gested by Mortonson and Hu [|| and we derive constraints 
on the neutrino mass in this different theoretical frame- 
work. It has been shown recently [6[ that a general reion- 
ization scenario can drastically alter the conclusions on 
inflationary parameters as the scalar spectral index n, 
putting its value in better agreement with the expecta- 
tions of a Harrison-Zel'dovich n = 1, spectrum. So 
it is definetely timely to investigate what kind of impact 
a general reionization scenario can have on the current 
CMB neutrino bound. The paper is organized as fol- 
lows: in the next Section we briefly describe the reion- 
ization parametrization assumed and the decomposition 
in principal components. We also describe the data anal- 
ysis method. In Section III we present the results of our 
analysis and in Section IV we discuss our conclusions. 



ANALYSIS METHOD 

We adopt the method, developed in Ref. [8|, based 
on principal components that provide a complete ba- 
sis for describing the effects of reionization on large- 
scale .E-mode polarization. Following Ref. j^, one can 
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parametrize the reionization history as a free function of 
redshift by decomposing x e {z) into its principal compo- 
nents: 

M 

where the principal components, Sfj,(z), are the eigen- 
functions of the Fisher matrix that describes the depen- 
dence of the polarization spectra on x e (z) (again, see 
Ref. Q), are the amplitudes of the principal com- 
ponents for a particular reionization history, and x[(z) 
is the WMAP fiducial model at which the Fisher matrix 
is computed and from which the principal components 
are obtained. In what follows we use the publicly avail- 
able S/j.(z) functions and varied the amplitudes for 
fi = 1,...,5 for the first five eigenfunctions. The eigen- 
functions are computed in 95 bins from redshift z min = 6 
to redshift z max — 30 with x e (z) = 1.08 for z < 3, 
x e {z) = 1.0 for 3 < z < 6 and x e (z) = 10~ 4 for z > 30.. 
Hereafter we refer to this method as the MH (Mortonson- 
Hu) case. 

We have then modified the Boltzmann CAMB code [9| 
incorporating the generalized MH reionization scenario 
as in [8j and extracted cosmological parameters from cur- 
rent data using a Monte Carlo Markov Chain (MCMC) 
analysis based on the publicly available MCMC package 
cosmomc flpj ] . 

We consider here a flat ACDM universe described by 
a set of cosmological parameters 

{w b ,w c ,w,,e s ,n,log[10 10 A s ]}, (2) 

where uji = fl^h 2 and w c = Q c h 2 are the physical baryon 
and cold dark matter densities relative to the critical den- 
sity, w v is the physical energy density in massive neutri- 
nos, S is the ratio of the sound horizon to the angular 
diameter distance at decoupling, A s is the amplitude of 
the primordial spectrum, and n is the scalar spectral in- 
dex. We assume 3 degenerate, massive neutrinos with 
the same mass: 

m v = 30.8eV" x oj v (3) 

In what follows we will use as standard parameter the 
value £m„ = 3m„. 

The extra parameters needed to describe the reioniza- 
tion are the five amplitudes of the eigenfunctions for the 
MH case and one single common parameter, the optical 
depth r, for the sudden reionization case. 

Our basic data set is the seven-year WMAP data [l| 
(temperature and polarization) with the routine for com- 
puting the likelihood supplied by the WMAP team. 

RESULTS 
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TABLE I. 95% c.l. errors on cosmological parameters in the 
case of sudden reionization and MH reionization. The upper 
limit on the neutrino mass is relaxed by ~ 43%. 



reionization scenario. As we can see from the table, the 
CMB constraint on the neutrino mass is weakened by 
~ 40% when a more general reionization scenario is con- 
sidered. This is not simply due to an increase in the 
parameter space but also due to degeneracies present be- 
tween the cosmological parameters. Considering the MH 
reionization scenario renders values of the spectral in- 
dex n in better agreement with the Harrison-Zel'dovich 
n = 1 value (see 0). This changes the relative am- 
plitude of the peaks in the CMB angular spectrum and 
makes models with higher neutrino mass more consistent 
with the WMAP data. Introducing a neutrino mass has 
indeed the effect of decreasing the gravitational poten- 
tial at recombination, increasing the small scale CMB 
anisotropjQ- This can be counterbalanced by decreasing 
the value of the spectral index n as clearly shown by the 
anti-correlation in the n-Y>m v plane. A general reioniza- 
tion scenario brings higher values of n in agreement with 
observations, immediately resulting in a better compati- 
bility of larger neutrino masses. It is worth noticing that 
while in the standard reionization scenario HZ spectra are 
excluded at about three standard deviations when mas- 
sive neutrinos are included in the analysis, in the MH 
case the n = 1 spectra are well consistent with the data 
and inside the la c.l. also with Sm^ ~ 0.5eV. 

In Figure [1] we show the constraints on the Sto^ vs 
n plane, while in Figure [2] we show the constraints on 
the Em^ vs erg plane. The filled contours assume MH 
reionization while the empty contours assume standard, 
sudden, reionization. As we can see, MH reionization 
allows for values of the spectral index n closer to 1 (as 
already pointed out in 6]), for a larger neutrino mass 
and for a lower amplitude. It is interesting to note 
that a neutrino mass can in principle accomodate lower 
values of erg with CMB data. When MH reionization 
is assumed even lower values of as are consistent with 



In Table I we compare the constraints on several cos- 
mological parameters in the case of standard or MH 



1 The effect of neutrino mass on CMB lensing for the WMAP data 
is negligible. 
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CMB data. A low value of as ~ 0.77 is preferred by the 
recent detection of diffuse Sunyaev-Zel'dovich effect by 
the South Pole Telescope [ll| experiment. 
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FIG. 1. Constraints on the Sm„ vs n plane. The filled con- 
tours assume MH reionization while the empty contours as- 
sume standard, sudden, reionization. 
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FIG. 3. Constraints on the Em„ vs Q m plane. The filled 
contours assume MH reionization while the empty contours 
assume standard, sudden, reionization. 
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FIG. 4. Constraints on the Em„ vs Q b h 2 plane. The filled 
contours assume MH reionization while the empty contours 
assume standard, sudden, reionization. 



FIG. 2. Constraints on the Em„ vs as plane. The filled 
contours assume MH reionization while the empty contours 
assume standard, sudden, reionization. 



Moreover, correlations exist with the matter density 
Q m , as we show in Figure [3] and (even if less pronounced) 
with the baryon physical density Q(,h 2 as we show in 
Figure |U 



CONCLUSIONS 

In conclusion, the details of the reionization processes 
in the late universe are not very well known. In the 



absence of a precise, full redshift evolution of the ion- 
ization fraction during the reionization period, a sim- 
ple parametrization, with a single parameter z r , has be- 
come the standard reionization scheme in numerical anal- 
yses. However, more general reionization scenarios are 
certainly plausible and their impact on the cosmological 
constraints should be carefully explored. 

In this Brief Report we have investigated the stabil- 
ity of the CMB constraints on neutrino masses in gen- 
eralized reionization scenarios. We have found that a 
more general treatment of reionization could potentially 
weaken the current CMB upper limit on £m„ by ~ 40%. 
Cosmological information from BAO for example can be 
added in order to reduce the uncertainty on the neutrino 
mass. However the lack of knowledge on dark energy 
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and the assumption made with regards to the equation 
of state could again affect the neutrino mass limit with 
large-scale structure data. Future data expected from the 
Planck [T3| satellite on large angular scale CMB polar- 
ization will help in clarifying the thermal history of the 
Universe and in ruling out exotic reionization scenarios 
that are still in agreement with present-day observations 
with WMAP. 
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